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Abstract 



We present an analysis of CMB anisotropies induced by global scalar fields 

^> ' in the large N limit. In this limit, the CMB anisotropy spectrum can be 

■ determined without cumbersome 3D simulations. We determine the source 
, functions and their unequal time correlation functions and show that they are 

CN I quite similar to the corresponding functions in the texture model. This leads 

' us to the conclusion that the large N limit provides a 'cheap approximation' 

Q> . to the texture model of structure formation. 

Ph: PACS numbers: 98.80. Cq 98.65.Dx 

2 ■ The anisotropies in the cosmic microwave background (CMB) have become an ex- 
c/3 . tremely valuable tool for cosmology. There are hopes that the measurements of the 
^ ' CMB anisotropy spectrum might lead to a determination of cosmological parame- 
ters like ^0, -f^O) ^B, A to within a few percent. The justification of this hope lies 
^ ' to a big part in the simplicity of the theoretical analysis. Fluctuations in the CMB 

■ can be determined almost fully within linear cosmological perturbation theory and 
are not severely influenced by nonlinear physics. 

Presently there are two competing classes of models which lead to a Harrison 
Zel'dovich spectrum of fluctuations: Perturbations may be induced during an infla- 
tionary epoch or they may be due to scaling seeds like, e.g., a self ordering global 
scalar field or cosmic strings (for a general definition of 'scaling seeds' see [Q). In 
the first class, the linear perturbation equations are homogeneous. In the second 
class they are inhomogeneous, with a source term due to the seed. The evolution of 
the seed is in general non-linear and complicated and therefore much less accurate 
predictions have been made so far for models where perturbations are induced by 
seeds. 

In this communication we discuss an especially simple model with seeds where 
the equation of motion for the seed perturbations can be solved explicitly. We 
consider a A^— component real scalar field (p with 0{N) symmetric potential V, 
which at T = is given by = A(0^ — 77^)^. At low temperatures, T ^ 77, can 
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be regarded as constrained to an — 1 sphere with radius r]. The scalar field then 
evolves according to the non-linear cr-model which is entirely scale free. In terms 
of the dimension- less variable [5 = (p/rj we find 

ap-{p.Dp)p = o , (1) 

with the condition (3"^ = 1. The non-linearity in this equation, —{f3 ■ Of3)(3 = 
{dfj,f5 ■ d^j3)l3 contains a sum over components. In the limit oo, this sum can 
be replaced by an ensemble average and the resulting linear equation of motion can 
be solved exactly. One obtains 

mt) = Atl:^f3Uk). (2) 

The index u is determined by the background matter model and varies between 
z/ = 2 in a radiation dominated background and i/ = 3 in a matter dominated 
background. The pre-factor A is chosen to ensure jS"^ = 1, 

' 15/16, ifz/ = 2, 
2835/128, ifz/ = 3. 

The components of Pin are assumed to be independent, Gaussian-distributed random 
variables with vanishing mean and dispersion {{f3in)'j) = 1/iV for all values of j. 
(Clearly, the variables {Pin)j cannot be completely independent, since they obey the 
condition Ei(An)f = 1-) 

Once the scalar field jS is known, we can calculate its energy-momentum tensor, 
the induced gravitational field and and its action on matter and radiation within 
linear cosmological perturbation theory. As has been discussed in 0, the energy 
density of a four component global scalar field is already quite close to the large 
A^ limit and there are thus justified hopes, that this simple model might provide 
a quite sensible approximation to the texture scenario for structure formation. On 
the other hand, we know that non-linearities, which lead to the mixing of scales 
and to the deviations from a Gaussian distribution, are crucial for some qualitative 
properties of defect models, like decoherence P, In the large— A^ limit, the only 
non-linearities are the quadratic expressions of the energy momentum tensor, and 
thus effects like decoherence might be mildened substantially in this model. 

This communication is an outline of a longer paper in preparation. We first 
briefly repeat the basic equations for the determination of the CMB anisotropy 
spectrum in the presence of seeds, and solve the equations (for given seed functions) 
in a simplified situation. We then outline the calculation of the relevant correlation 
functions and present some results. We end with conclusions and an outlook. 

The coefficients of the angular power spectrum of CMB anisotropies are related 
to the 2-point function according to 
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For pure scalar perturbations, neglecting Silk damping, the C^'s are given by[]I[ 
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A3,(A;, t^ccMkto) - Vr{k, tdec)j;(Ho) + A; / - ^){k, t')f,{k{to - t'))dt'. 

(5) 

For large £ this spectrum is corrected by Silk damping, which can be approxi- 
mated by multiplying /S.f (k) with an exponential damping envelope 0]. 

We want to consider the situation where fluctuations are induced by seeds. We 
restrict ourselves to scalar perturbations. The energy momentum tensor of scalar 
seed perturbations can be parameterized by the following four functions: fp the 
energy density of the seed, fp the pressure of the seed, a potential for the energy 
flux of the seed and the potential for anisotropic stresses of the seed [^, 0] (see 
below). The linear cosmological perturbation equations are then of the form 
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where P is a first order linear differential operator, X is a vector consisting of all, say 
m, gauge invariant perturbation variables of the cosmic fluids (like A^^, and, e.g. 
the corresponding variables for the cold dark matter (CDM), ...). F = {fp, fp, f^, /jr) 
is the source vector and M is a, in general time dependent, m x 4 matrix. 
The general solution to Eq. (0) is of the form 



G){t,t')F,{t')dt' 
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where G is the Green's function of the differential operator V. The Bardeen poten- 
tials $ and \l/ are algebraic combinations of the fluid variables Xj and the source 
functions Fj, 

- m) = pKt)x,{t) + Q\m{t) . (8) 

Inserting this solution in Eq. (13) we obtain 
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The expectation value in Eq. (^ thus consists of time integrals of unequal time 
correlations of the source functions, which are of the generic form 
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To calculate the C^'s we need therefore need to know the unequal time correlations of 
the seed functions Fi and the Green's functions for the cosmological model specified. 
In general this is a quite formidable task. Here we shall just discuss a toy model. A 
somewhat more complicated example is given in 0. 

We consider a pure radiation universe with vanishing spatial curvature. In this 
case, the linear perturbation equations are given by 
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where x = kt and a prime denotes a derivative w.r.t x. The energy momentum 
tensor of the source enters via the combinations and $5 = e^fp/k"^ + 3fv/{kx)). 
These equations can be combined to a second order differential-equation for A^^ 
alone, 
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leading to the Greens function 
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On super-horizon scales (x -C 1) the solutions of the homogeneous equations consist 
of one constant and one decaying, oc 1/x, mode, while for x ^ 1 we obtain two 



4 



oscillating modes. The general solution with source term S{x) and initial condition 
A^^(O) = Vr{0) = is 

Agr{x) = r G\x,x')S{x')dx' (17) 
Jo 

Vr{x) = r G^{x,x')S{x')dx' where (18) 
Jo 

G' = G and G' = . (19) 

4 ax 

Together with Eq. ([TI|) and (|12D it is now straight forward to determine the integral 
kernels A,B,C,D,E and H and thus, for given source correlation functions, the C^'s. 

Let us therefore discuss the source correlation functions of scalar field sources in 
the large limit. The seed functions are given by 0]: 
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Using the fact that the initial fields are uncorrelated and Gaussian distributed 
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and using the exact solution Eq. (H) we find the power spectra and the unequal time 
correlation functions of the seed variables Fi. Below we give explicit expressions for 
the power spectra of fp , and and, as an example, the unequal time correlation 
function for f^j . These integrals can be evaluated numerically, examples of which are 
shown in Figs. 1 and 2. We define x = tk and we set 
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where we have set r = t' /t in the last equation. The behavior of these functions 
on very large and very small scales can be obtained analytically. On super hori- 
zon scales, X — i> 0, the power spectra for fp , fp and behave like white noise. 
Numerically we have found 
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From general arguments P, IH, we would have expected also to behave like 
white noise on super-horizon scales. However, from Eq. ([28| ) we find that f^^ diverges 
at small x like l/x"^. Even though we do not quite understand this result, it does 
not lead to divergent Bardeen potentials, if we allow for anisotropic stresses in the 
matter (like e.g. from a component of massless neutrinos). In this case it can be 
shown [0] that compensation arranges the anisotropic stresses in the fluid, pU, such 
that + pU oc x'^Jt^. Therefore, also the anisotropic stresses contribute a white 
noise component to the Bardeen potentials on super-horizon scales, namely: 
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In the limit x 



oo the source functions decay like 
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In Fig. 1 we plot t < |/pp > (x) as obtained from Eq. (|26|) , and compare it to the 
corresponding function found by 3D simulations of the texture model. 
The normalized unequal time correlation functions are defined by 



In the large N limit, the correlation functions decay like power laws. For r = t'/t 
we find in the limit r ^ 1 , kt' ^ 1 the behavior Ci oc r~'^', with 



It is not quite clear to us, whether this behavior is reproduced in the texture model. 
Due to the arguments given at the beginning, it may well be that = 4 and N ^ oo 
show a different decoherence behavior. Originally (taking into account the numerical 
accuracy of about 10% of the 3D simulations) we approximated the decoherence in 
the texture model with an exponential decay law. However, comparing the unequal 
time correlation functions for shown in Figs. 2 and 3 for the large limit and a 3D 
simulation of the texture model respectively, we realize, that they agree extremely 
well in the numerically most reliable, central region, and the seemingly stochastic 
higher order oscillations also found in the texture model might actually be real (see 
Fig. 4), leading to power law decoherence. 

Using these source functions, we have determined the CMB anisotropy spectrum 
induced by the large N limit of a self ordering scalar field for a spatially flat cosmo- 
logical model with CDM, radiation and baryons. Since decoherence is so weak for 
large A^, we used the approximation of perfect coherence, Ci = 1. This simplification 
has been used so far for all 'analytic approximations' of 0{N) models and, e.g., in 
the case of textures, = 4, it seems to agree reasonably with numerical simula- 
tions This will certainly be even more so in the large A^ limit. The influence of 
decoherence on the CMB power spectrum is discussed in [|]. 

We have shown that the large A^ limit of global scalar fields provides a model 
of seeded structure formation where CMB anisotropies can be determined without 
cumbersome numerical simulations and thus with much higher accuracy and large 
dynamical range at relatively modest costs. Determining the correlation functions of 
the seed variables just requires numerical convolution of Bessel functions multiplied 
with powers. Once the seed correlation functions are known, perturbations in matter 
and radiation can be calculated by solving a system of linear perturbation equations, 
very similar to the homogeneous case of inflationary perturbations. 
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7p = 3/2, 7p = 3/2, 7. = 3/2, 7. = 5/2 . 



(40) 



We believe that the large limit has many features in common with the texture 
model of structure formation and thus provides a "cheap approximation" to this 
model. The most obvious difference between the analytic limit and the texture 
model is the decoherence behavior. In the large N limit, the field evolution is 
linear and non-linearities, which are responsible for decoherence, enter only via the 
quadratic expressions of the energy momentum tensor. We thus expect decoherence 
to be somewhat weaker in the large N limit. 

In a forthcoming paper, we plan to work out the large limit in more detail, and 
to study the dependence of the resulting CMB anisotropy spectrum of cosmological 
parameters. We also want to investigate more fully the comparison of the large N 
source functions with the source functions found in 3D simulations of the texture 
model. The limit discussed here provides an very useful toy model for structure 
formation with scaling seeds for which decoherence is not important. 
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Figure 1: The source function fp-. In the large N hmit (full line), the approximation 
used to calculate the C^s in Fig. 5 (dashed line) and from a 3 dimensional numerical 
computation (dot-dashed line) for the texture model (A^ = 4). 
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Figure 2: The unequal time correlation function for = fp + Sfp at fixed t as 
function of t' and k in the large limit. Negative values are set to zero. 
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Figure 3: The same as Fig. 2 for the texture model. The similarity is obvious. The 
high order oscillations which are very pronounced in the large limit are washed 
out or absent in the texture model. Whether this is a real feature or just numerical 
inaccuracy or both is not yet clear. 




Figure 4: A cut through Fig. 2 (solid line) and Fig. 3 (dashed line) at kt = 3.9. The 
central peaks are in very good agreement. Secondary peaks do not agree and the 
decay law for the texture model is difficult to predict from this data. 
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Figure 5: The CMB anisotropy spectrum, obtained by using polynomial fits for 
the source functions as shown in Fig. 1. Only scalar perturbations are included. 
The Sachs Wolfe part is indicated by the dotted line, the dashed line represents the 
acoustic contributions. Silk damping is not included. 
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